Coexistence of attractors in a laser diode with optical feedback from a large external cavity

The coexistence of several attractors in the coherent collapsed state of a laser diode with optical feedback is investigated. These attractors are unstable tori that bifurcate from different external cavity modes. As the feedback rate is increased, these tori undergo different types of quasiperiodic routes, such as frequency locking, period doubling, or the appearance of a third incommensurate frequency. In the fully developed coherent collapsed state, the tori are all unstable, and the phenomenon of intermittence appears. Lyapunov exponent calculations demonstrate that this state presents hyperchaotic and high-dimensional dynamics. These results are explained qualitatively in terms of the multiattractor behavior found.Peer ReviewedPostprint (published version

Distribution of residence times of time-delayed bistable systems driven by noise

I study bistable time-delayed feedback systems driven by noise. Based on a two-state model with transition rates depending on the earlier state of the system I calculate analytically the residence-time distribution function. I show that the distribution function has a detailed structure, reflective of the effect of the feedback. By using an adequate indicator I give evidence of resonant behavior in dependence on the noise level. I also predict that this feedback-induced effect might be observed in two well-known optical bistable systems.Peer ReviewedPostprint (published version

Outlier Mining Methods Based on Graph Structure Analysis

Outlier detection in high-dimensional datasets is a fundamental and challenging problem across disciplines that has also practical implications, as removing outliers from the training set improves the performance of machine learning algorithms. While many outlier mining algorithms have been proposed in the literature, they tend to be valid or efficient for specific types of datasets (time series, images, videos, etc.). Here we propose two methods that can be applied to generic datasets, as long as there is a meaningful measure of distance between pairs of elements of the dataset. Both methods start by defining a graph, where the nodes are the elements of the dataset, and the links have associated weights that are the distances between the nodes. Then, the first method assigns an outlier score based on the percolation (i.e., the fragmentation) of the graph. The second method uses the popular IsoMap non-linear dimensionality reduction algorithm, and assigns an outlier score by comparing the geodesic distances with the distances in the reduced space. We test these algorithms on real and synthetic datasets and show that they either outperform, or perform on par with other popular outlier detection methods. A main advantage of the percolation method is that is parameter free and therefore, it does not require any training; on the other hand, the IsoMap method has two integer number parameters, and when they are appropriately selected, the method performs similar to or better than all the other methods tested.Peer ReviewedPostprint (published version

Correlation lags give early warning signals of approaching bifurcations

Identifying approaching bifurcations and regime transitions from observations is an important challenge in time series analysis with practical applications in many fields of science. Well-known indicators are the increase in spatial and temporal correlations. However, the performance of these indicators depends on the system under study and on the type of approaching bifurcation, and no indicator provides a reliable warning for any system and bifurcation. Here we propose an indicator that simultaneously takes into account information about spatial and temporal correlations. By performing a bivariate correlation analysis of signals recorded in pairs of adjacent spatial points, and analyzing the distribution of lag times that maximize the cross-correlation, we find that the variance of the lag distribution displays an extreme value that is a consistent early warning indicator of the approaching bifurcation. We demonstrate the reliability of this indicator using different types of models that present different types of bifurcations, including local bifurcations (transcritical, saddle-node, supercritical and subcritical Hopf), and global bifurcations.This work was funded by the Spanish Ministerio de Ciencia, Innovación y Universidades (PGC2018-099443-B-I00 ) and the ICREA ACADEMIA program of Generalitat de Catalunya.Peer ReviewedPostprint (published version

Entropy-based early detection of critical transitions in spatial vegetation fields

In semiarid regions, vegetated ecosystems can display abrupt and unexpected changes, i.e., transitions to different states, due to drifting or time-varying parameters, with severe consequences for the ecosystem and the communities depending on it. Despite intensive research, the early identification of an approaching critical point from observations is still an open challenge. Many data analysis techniques have been proposed, but their performance depends on the system and on the characteristics of the observed data (the resolution, the level of noise, the existence of unobserved variables, etc.). Here, we propose an entropy-based approach to identify an upcoming transition in spatiotemporal data. We apply this approach to observational vegetation data and simulations from two models of vegetation dynamics to infer the arrival of an abrupt shift to an arid state. We show that the permutation entropy (PE) computed from the probabilities of two-dimensional ordinal patterns may provide an early warning indicator of an approaching tipping point, as it may display a maximum (or minimum) before decreasing (or increasing) as the transition approaches. Like other spatial early warning indicators, the spatial permutation entropy does not need a time series of the system dynamics, and it is suited for spatially extended systems evolving on long time scales, like vegetation plots. We quantify its performance and show that, depending on the system and data, the performance can be better, similar or worse than the spatial correlation. Hence, we propose the spatial PE as an additional indicator to try to anticipate regime shifts in vegetated ecosystems.Peer ReviewedPostprint (published version

Sub-threshold signal encoding in coupled FitzHugh-Nagumo neurons

Despite intensive research, the mechanisms underlying the neural code remain poorly understood. Recent work has focused on the response of a single neuron to a weak, sub-threshold periodic signal. By simulating the stochastic FitzHugh-Nagumo (FHN) model and then using a symbolic method to analyze the firing activity, preferred and infrequent spike patterns (defined by the relative timing of the spikes) were detected, whose probabilities encode information about the signal. As not individual neurons but neuronal populations are responsible for sensory coding and information transfer, a relevant question is how a second neuron, which does not perceive the signal, affects the detection and the encoding of the signal, done by the first neuron. Through simulations of two stochastic FHN neurons we show that the encoding of a sub-threshold signal in symbolic spike patterns is a plausible mechanism. The neuron that perceives the signal fires a spike train that, despite having an almost random temporal structure, has preferred and infrequent patterns which carry information about the signal. Our findings could be relevant for sensory systems composed by two noisy neurons, when only one detects a weak external input.Peer ReviewedPostprint (published version

Neuronal coupling benefits the encoding of weak periodic signals in symbolic spike patterns

The biophysical mechanisms by which an input signal elicits a neuronal response are well known (sufficiently large inputs change the membrane potential of the neuron and generate electrical pulses, known as action potentials or spikes), yet, a good understanding of how neurons use these spikes to encode the signal information remains elusive. Recent theoretical studies have focused on how neurons encode a weak periodic signal (that by itself is unable to generate spikes) in a noisy environment, where stochastic electrical fluctuations that do not encode any information occur. Analyzing spike sequences generated by individual neurons and by two coupled neurons (that were simulated with the stochastic FitzHugh–Nagumo model), it has been found that the relative timing of the spikes can encode the signal information. Using a symbolic method to analyze the spike sequence, preferred and infrequent spike patterns were detected, whose probabilities vary with both, the amplitude and the frequency of the signal. To investigate if this encoding mechanism is plausible also for neuronal ensembles, here we analyze the activity of a group of neurons, when they all perceive a weak periodic signal. We find that, as in the case of one or two coupled neurons, the probabilities of the spike patterns, now computed from the spike sequences of all the neurons, depend on the signal’s amplitude and period, and thus, the patterns’ probabilities encode the information of the signal. We also find that the resonances with the period of the signal or with the noise level are more pronounced when a group of neurons perceive the signal, in comparison with when only one or two coupled neurons perceive it. Neuronal coupling is beneficial for signal encoding as a group of neurons is able to encode a small-amplitude signal, which could not be encoded when it is perceived by just one or two coupled neurons. Interestingly, we find that for a group of neurons, just a few connections with one another can significantly improve the encoding of small-amplitude signals. Our findings indicate that information encoding in preferred and infrequent spike patterns is a plausible mechanism that can be employed by neuronal populations to encode weak periodic inputs, exploiting the presence of neural noise.Peer ReviewedPostprint (author's final draft

Dynamics of a semiconductor laser with feedback and modulation: experiments and model comparison

© 2222 Optica Publishing Group. Users may use, reuse, and build upon the article, or use the article for text or data mining, so long as such uses are for non-commercial purposes and appropriate attribution is maintained. All other rights are reserved.We study experimentally and numerically the dynamics of a semiconductor laser near threshold, subject to optical feedback and sinusoidal current modulation. The laser operates in the low frequency fluctuation (LFF) regime where, without modulation, the intensity shows sudden spikes at irregular times. Under particular modulation conditions the spikes lock to the modulation and their timing becomes highly regular. While the modulated LFF dynamics has received a lot of attention, an in-depth comparison with the predictions of the Lang-Kobayashi (LK) model has not yet been performed. Here we use the LK model to simulate the laser dynamics and use the Fano factor to quantify the regularity of the timing of the spikes. The Fano factor is calculated by counting the number of spikes in successive segments of the intensity time-series and keeps information about temporal order in the spike sequence that is lost when the analysis is based on the distribution of inter-spike intervals. Here we compare the spike timing regularity in experimental and in simulated spike sequences as a function of the modulation amplitude and frequency and find a good qualitative agreement. We find that in both experiments and simulation for appropriate conditions the spike timing can be highly regular, as revealed by very small values of the Fano factor.Institució Catalana de Recerca i Estudis Avançats (Academia); Ministerio de Ciencia, Innovación y Universidades (PGC2018-099443-B-I00).Peer ReviewedPostprint (published version

Quantifying the synchronization of the spikes emitted by coupled lasers

Synchronization phenomena is ubiquitous in nature, and in spite of having been studied for decades, it still attracts a lot of attention as is still challenging to detect and quantify, directly from the analysis of noisy signals. Semiconductor lasers are ideal for performing experiments because they are stochastic, nonlinear, and inexpensive and display different synchronization regimes that can be controlled by tuning the lasers’ parameters. Here, we analyze experiments done with two mutually optically coupled lasers. Due to the delay in the coupling (due to the finite time the light takes to travel between the lasers), the lasers synchronize with a lag: the intensity time traces show well-defined spikes, and a spike in the intensity of one laser may occur shortly before (or shortly after) a spike in the intensity of the other laser. Measures that quantify the degree of synchronization of the lasers from the analysis of the intensity signals do not fully quantify the synchronicity of the spikes because they also take into account the synchronization of fast irregular fluctuations that occur between spikes. By analyzing only the coincidence of the spike times, we show that event synchronization measures quantify spike synchronization remarkably well. We show that these measures allow us to quantify the degree of synchronization and, also, to identify the leading laser and the lagging one.C.M. acknowledges the support of Ministerio de Ciencia, Innovación y Universidades (No. PID2021-123994NB-C21) and Institució Catalana de Recerca i Estudis Avançats (Academia).Peer ReviewedPostprint (author's final draft

Polarization dynamics in vertical-cavity surface-emitting lasers with optical feedback through a quarter-wave plate

Square-wave switching of the intensities of the orthogonal linearly polarized components of the output of a vertical-cavity surface-emitting laser(VCSEL) found previously in experiments with polarization-changing optical feedback, is not found in rate equation models incorporating only birefringence and gain anisotropy, but is found in the model for VCSELs developed by San Miguel, Feng, and Moloney [M. San Miguel, Q. Feng, and J. V. Moloney, Phys. Rev. A 52, 1729 (1995)]. The dynamics is sensitive to both the feedback strength and the relaxation rate of the magnetization in the quantum well sublevels.Peer ReviewedPostprint (published version